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ABSTRACT 

O 

o 

CN ' We report the detection of differential gas column densities in three gravita- 

tional lenses, MG 04f4+0534, HE 1104-1805, and PKS 1830-211. Combined 
with the previous differential column density measurements in B 1600+434 and 
Q 2237+0305 and the differential extinction measurements of these lenses, we 
probe the dust-to-gas ratio of a small sample of cosmologically distant normal 
galaxies. We obtain an average dust-to-gas ratio of E(B — V)/ N^= (1.4 ± 0.5) x 
10 -22 mag cm 2 atoms -1 with an estimated intrinsic dispersion in the ratio of 



> 



~ 40%. This average dust-to-gas ratio is consistent with the average Galactic 
t^. ■ value of 1.7 x 10 -22 mag cm 2 atoms -1 and the estimated intrinsic dispersion is 



also consistent with the 30% observed in the Galaxy. 
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i-£h ■ Subject headings: dust-to-gas ratio 
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^ ' 1. Introduction 



Gravitational lenses are valuable tools in many astrophysical applications. One of them 
is to study the properties of the interstellar medium (ISM) in the lens galaxies at the locations 
where we observe the quasar images. By studying the differences in the absorption (or 
extinction) between multiple images, any contamination from the Galaxy in the foreground 
or the quasar host galaxy in the background is essentially eliminated and we can be confident 
that we are probing the ISM of the lens galaxy. Several groups (e.g., Nadeau et al. 1991; 
Falco et al. 1999; Toft, Hjorth & Burud 2000; Motta et al. 2002; Wucknitz et al. 2003; 
Munoz et al. 2004; Mediavilla et al. 2005) have studied dust extinction and the dust 
extinction law in lens galaxies. In X-rays, Dai et al. (2003) and Dai & Kochanek (2005) 
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measured the differential column densities in two lenses Q 2237+0305 and B 1600+434 using 
a similar method. By combining the X-ray and optical measurements of the differences in 
gas and dust properties between the lensed images, we can estimate the dust-to-gas ratio 
AE(B — V)/ANn of the lens galaxies assuming the differences between the extinction and 
gas absorption are due to the same parcel of the ISM. 

The dust-to-gas ratio, normally calculated from E(B — V)/Nn, is a basic property of 
the ISM, but our direct knowledge of dust-to-gas ratios is largely based on our own Galaxy. 
In particular, Bohlin, Savage & Drake (1978) measured the dust-to-gas ratios of the ISM 
towards 100 stars in the Milky Way and found a correlation of EiB — V)/Nu = 1.7 x 10~ 22 
mag cm 2 atoms -1 with typical scatter of ~30%. Subsequent studies have extended these 
measurements to additional lines of sight in our Galaxy and to the LMC and SMC, but 
the basic results are little altered aside from some evidence that the dust-to-gas ratios of 
the LMC and SMC are modestly lower (~ 30%) than in the Galaxy (see the review by 
Draine 2003 and references therein). However, because of the difficulties in determining the 
absorption column densities and the extinction of distant galaxies, little is known about the 
dust-to-gas ratios of their ISM. Just as with studies of extinction laws, gravitational lenses 
provide the first potential probe for studying the dust-to-gas ratios of cosmologically distant 
galaxies (0 < z < 1) with the precision of studies of the ISM in the Galaxy. Here we measure 
the differential column density of three additional lenses and combine them with existing 
measurements for two lenses to estimate the average dust-to-gas ratio of the lens galaxies. 



2. Sample Selection and Data Reduction 

Over twenty gravitational lenses have been observed in X-rays with the Chandra X- 
ray Observatory (Weisskopf et al. 2002). From these systems, we selected the seven sys- 
tems from the Falco et al. (1999) survey of extinction in lens galaxies with differential 
extinction measurements exceeding AE(B — V) > 0.05 mag. For two of these systems, 
LBQS 1009—0252 and MG 2016+112, the Chandra spectra of the individual images contain 
too few photons to allow an analysis of differential absorption, and the differential absorp- 
tions in B 1600+434 and Q 2237+0305 have been reported previously (Dai et al. 2003; 
Dai & Kochanek 2005). Here we analyze the remaining three lenses, MG 0414+0534, 
HE 1104-1805, and PKS 1830-211. 

MG 0414+0534 and HE 1104-1805 were observed with the Advanced CCD Imaging 
Spectrometer (ACIS, Garmire et al. 2003) on board Chandra for 97 and 49 ks on January 
8, 2002 and June 10, 2000, respectively. PKS 1830—211 was observed with the High Energy 
Transmission Grating Spectrometer (HETGS, Canizares et al. 2005) and ACIS on board 
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Chandra for 51 ks on June 25, 2001. The Chandra data were reduced with the CIA0 software 
tools provided by the Chandra X-ray Center (CXC) following the standard threads on the 
CXC website. 1 Only events with standard ASCA grades of 0, 2, 3, 4, and 6 were used in the 
analysis. We improved the image quality of the data by removing the pixel randomization 
applied to the event positions by the standard pipeline. In addition, we applied a subpixel 
resolution technique (Tsunemi et al. 2001; Mori et al. 2001) to the events on the S3 chip 
of ACIS where the lensed images are located. 



We extracted spectra for each of the lensed quasar images. Although the PKS 1830—211 
observation was a grating observation, we simply analyzed the spectra obtained from the 
zeroth order image to extract and study spectra of the two quasar images separately. We did 
not analyze the grating spectrum of PKS 1830—211 because it cannot resolve the spectra of 
the two images, which is essential in this study. We fitted the spectra of the lensed quasars 
using XSPEC VI 1 .3. 1 (Arnaud 1996) over the 0.3-8 keV observed energy range. The X-ray 
spectrum of the ith image can be generally modeled as 



where Ni(E, t) is the number of photons per unit energy interval, iVn gwj A% j, and Nn^ Src are 
the equivalent hydrogen column density at our Galaxy, the lens galaxy at image i, and the 
source galaxy, cr(E) is the photo-electric absorption cross-section, and Atj is the time-delay. 
The difference between the absorption at the lens, AA^h, can be obtained by analyzing the 
individual spectra of the lensed images. We neglected the time-delay between the images 
and the absorption at the source redshift in this analysis. We used the standard wabs and 
zwabs models in XSPEC to model the Galactic absorption and absorption at the lens redshift. 
The wabs and zwabs models use cross sections from Morrison & McCammon (1983) and 
assume a solar elemental abundance from Anders & Ebihara (1982). 



3. Spectral Analysis 




exp{-o-(E)iV ] 



H,Ga/ - ° [E(l + Zl )\ A H , - a [E(l + z s )\ N^ Src {t - At,)} ,(1) 



1 The CXC website is at http://cxc.harvard.edu/. 
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3.1. MG 0414+0534 

MG 0414+0534, discovered by Hewitt et al. (1992), is a four-image lens system with 
a source redshift of z s = 2.639 (Lawrence et al. 1995) and a lens redshift of zi = 0.9584 
(Tonry & Kochanek 1999). The Chandra image of MG 0414+0534 was presented in Chartas 
et al. (2002). We extracted the spectra of the relatively well separated images A (Al and 
A2 combined), B and C from the 97 ks Chandra observation. The Al and A2 images are 
too close (0'.'4 apart) to obtain the individual spectra of the images. 

We modeled the spectra with a power-law modified by absorption at the lens and Galac- 
tic absorption. The Galactic absorption was fixed as A" H = 0.11 x 10 22 cm" 2 (Dickey & 
Lockman 1990). We fitted the three spectra simultaneously using the same power-law pho- 
ton index for all spectra but differing amounts of absorption at the lens redshift for each 
image. The results are presented in Table 1. Our estimate of the power-law photon index, 
T = 1.69^0 oi) i s consistent with the range of values from Chartas et al. (2002). Figure 1 
shows the spectra of the three images and their best fit models. It is clear from the spectra 
that the spectrum of image A is more heavily absorbed than that of images B and C. We 
estimate that the differences in the column densities are ATVha.b — (0.33 + 0.10) x 10 22 cm" 2 , 
AiV HA)C = (0.20 ± 0.15) x 10 22 cm" 2 and A# Hc ,b = 0.1^ x !0 22 cm" 2 . We used XSPEC to 
explore the full parameter space of the column densities in order to correctly estimate the 
uncertainties in the differential column densities. We also verified that adding absorption 
in the source does not alter the estimates of the differential column densities between the 
images in the lens galaxy. 

Falco et al. (1999) measured the differential extinction of the four images as AE(B — 
V) = 0.09, 0.31, 0.02 and mag for images Al, A2, B and C, respectively. The differential 
extinction between images B and C is very small. The differential extinction between image 



Table 1. Spectral Fitting Results 



Quasar 


ModeP 


Galactic JVh 
(10 22 cm- 2 ) 


Nn (Image A) 
(10 22 cm- 2 ) 


JVh (Image B) 
(10 22 cm- 2 ) 


Nn (Image C) 
(10 22 cm- 2 ) 


r 


Xlidof) 


MG 0414+0534 
HE 1104-1805 
PKS 1830-211 


wabs (zwabs (pow) ) 
wabs (zwabs (pow) ) 
wabs (zwabs (pow) ) 


0.11 (fixed) 

< 0.004 
0.22 (fixed) 


0.86 ±0.06 
0.06 ± 0.03 
1 0+ ' 4 


0.53 ±0.10 
< 0.01 
2.8 ±0.6 


64+ ' 22 


1.68 ±0.06 
1.02 ±0.06 


0.95(318) 
0.88(84) 
1.08(237) 



Note. — Simultaneous fits to the individual spectra of the lensed quasars. The spectra are constrained to have the same intrinsic 
power-law photon index but can have different absorption column densities at the redshift of the lens. 



a Thc wabs, zwabs, and pow represent the Galactic absorption, absorption at the lens, and the intrinsic power-law spectrum, respectively. 
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Al and A2 is large with AE(B — V) — 0.22 mag. We used the two brightest images to 
estimate the dust-to-gas ratio. Because we only analyzed the X-ray spectrum of the image 
A (Al and A2 combined), we used an average differential extinction of AE(B — V) — 
0.18 mag between images A and B in our dust-to-gas ratio estimate. Combined with our 
differential absorption measurement, we obtained a dust-to-gas ratio of E(B — V)/Nr = 
(0.55±0.33) x 10~ 22 mag cm 2 atoms -1 . The error-bar was estimated by replacing the average 
differential extinction with values for the Al and A2 images. This ambiguity dominates any 
interpretation of the differential absorption measurement. 



3.2. HE 1104-1805 



HE 1104—1805 is a two-image lens discovered by Wisotzki et al. (1993) with a source 
redshift of z s = 2.319 (Wisotzki et al. 1993) and a lens redshift of z\ = 0.729 (Lidman 
et al. 2000). The two images are well separated (AO = 3'.'2), so we can easily extract the 
individual spectra of both images. The combined spectrum of both images of HE 1104—1805 
was previously analyzed by Dai et al. (2004). In fitting the combined spectrum, as done by 
Dai et al. (2004), the low total absorption obtained indicates that adding a small component 
of Galactic absorption (Ah = 0.05 x 10 22 cm -2 ) suffices to fit the data with a power law index 
of r = 1.86to'o4 for the intrinsic spectrum. We again see the advantages of using differential 
measurements to eliminate ambiguities as to the location of the absorption. We obtained a 
differential column density between the two images of ANua,b = (0.055 ± 0.030) x 10 22 cm~ 2 
(Table 1) by placing the absorption at the lens redshift. Figure 2 shows the spectra of the 
two images and their best fit models. The power-law photon index, T = 1.68 ±0.06, obtained 
in this model is slightly harder than the values from Dai et al. (2004). Curiously, we find 
that image A shows more absorption than image B, which is the reverse of the differential 
extinction estimated by Falco et al. (1999). This is probably due to chromatic microlensing 
of the images by the stars in the lens galaxy being misinterpreted as extinction, but it 
could also be explained by giving the ISM very different dust-to-gas ratios at the positions 
of the two images. Recent observations of HE 1104—1805 (Wisotzki et al. 1993, 1995; 
Gomez- Alvarez et al. 2004) indicates that the chromaticity disappears when the microlensing 
perturbations decrease, supporting the chromatic microlensing interpretation. In addition, 
it is possible that the amount of ISM varied between the optical and X-ray observations. 
However, considering the short time-scale, it is unlikely that the ISM of a normal galaxy 
would vary significantly. In any case, we will be unable to include this lens in an estimate 
of the mean dust-to-gas ratio. 
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Fig. 1.— Spectra of images A (top), B and C (bottom) of MG 0414+0534. 
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Fig. 2.— Spectra of images A (top) and B (bottom) of HE 1104-1805. 
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3.3. PKS 1830-211 

PKS 1830—211 (Pramesh Rao & Subrahmanyan 1988; Subrahmanyan et al. 1990; 
Jauncey et al. 1991) consists of two z s = 2.507 (Lidman et al. 1999) quasar images separated 
by 1" and lensed by a z t = 0.886 (Wiklind & Combes 1996; Gerin et al. 1997) spiral galaxy 
(Winn et al. 2002). Using our standard method, but applying it to the zeroth order HETGS 
image, we find that ANub,a = 1.8t°;g x 10 22 cm -2 (Table 1) with the Galactic absorption 
fixed at iVn = 0.22 x 10 22 cm" 2 (Dickey & Lockman 1990). Figure 3 shows the spectra of 
the two images and their best fit models. We obtained an extremely hard power-law photon 
index, T = 1.02 ±0.06, considerably harder than the values obtained from most other radio- 
loud quasars ((r) ~ 1.5-1.6, e.g., Cappi et al. 1997; Reeves & Turner 2000) or blazars 
((r) ~ 2.2, e.g., Donato et al. 2001). However, it is consistent with the previous ROSAT 
result (r = 1+ 2 , Mathur & Nair 1997). 

We also estimated the absorption from the combined spectrum of the two images to 
compare with previous measurements. Again, we fixed the Galactic absorption at Nr = 
0.22 x 10 22 cm -2 . The best fit to the combined spectrum of Nu = (1-8 ± 0.4) x 10 22 cm -2 
coincidentally matches the differential absorption measurement. The combined absorption 
and the power-law photon index obtained, T = 1.01 ± 0.06, also agree with the analysis of 
the combined spectrum by de Rosa et al. (2005). This combined absorption is somewhat 
lower than the estimates from molecular absorption in the radio (N(H 2 ) ~ 2.5 x 10 22 cm" 2 , 
Wiklind k Combes 1996, 1998, Gerin et al. 1997) or the column of 3.5 x 10 22 cm" 2 required 
to match combined ROSAT spectrum of the two images with an assumed X-ray photon 
index of V = 2.2 (Mathur & Nair 1997). However, it is consistent with their measurement of 
Nu = l-7l 2 '5 x 10 22 cm ~ 2 allowing for a harder power-law photon index of T = it 2 (Mathur 
& Nair 1997). The very large absorption difference, ANn BA = 1.8* x 10 22 cm -2 , is 
consistent with the large differential extinction of AE(B — V) — 3.00 mag (Falco et al. 
1999) between images B and A. Combining these measurements, the dust-to-gas ratio is 
E(B — V)/N H = (1.7±0.6) x 10~ 22 mag cm 2 atoms -1 , consistent with the standard Galactic 
value. 



4. Discussion 

Table 2 summarizes the five systems for which we have estimates of both the differential 
extinction and the absorption column density. This includes the three lenses analyzed here 
(MG 0414+0534, HE 1104-1805 and PKS 1830-211) and the previously analyzed systems 
Q 2237+0305 (Dai et al. 2003) and B 1600+434 (Dai & Kochanek 2005). We did not 
include B 0218+357, as the dust-to-gas ratio estimated by Falco et al. (1999) was based on 
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Fig. 3.— Spectra of images A (top) and B (bottom) of PKS 1830-211. 



Table 2. The Dust-To-Gas Ratio of High Redshift (z > 0) Galaxies 











AJV H 


AE(B - V) 


E(B - V)/N H 


Lens 




type 


Between Imag 


cs (10 22 cm- 2 ) 


(mag) 


(10~ 22 mag cm 2 atoms -1 ) 



MG 0414+0534 


0.9584 


elliptical 


A, B 


0.33 ± 0.10 


0.18 a 


0.55 + 0.33 


HE 1104-1805 


0.729 


elliptical 


A, B 


0.055 ± 0.030 


-0.07 + 0.01 


b 


B 1600+434 


0.41 


spiral 


B, A 


26+ 17 


0.10 + 0.03 


0.38 + 0.25 c 


PKS 1830-211 


0.886 


spiral 


B, A 


1 s+ - 5 


3.00 + 0.13 


1.7 + 0.6 


Q 2237+0305 


0.0395 


spiral d 


A, C 


0.04 ± 0.03 


0.11 + 0.03 


2.8 + 2.2 



a The average extinction value of images Al and A2 was used in the calculation and the error-bar is estimated by replacing 
the average differential extinction with values for the Al and A2 images. 

b The dust-to-gas ratio was not estimated for HE 1104—1805 because the the differential extinction and TVjj column 
measurements give opposite signs. 

c The number provided in Dai & Kochanck (2005) is the gas-to-dust ratio rather than the dust-to-gas ratio. 

d Q 2237+0305 is lensed by the central bulge of a spiral galaxy. 
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Fig. 4. — AiVn versus AE(B — V) for the gravitational lenses (squares) and A^ H versus 
E(B — V) for Galactic stars from Bohlin et al. (1978, filled circles). The dashed line shows 
the mean Galactic relation. 
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very uncertain estimates (N(H 2 ) ~ (0.2 — 5.0) x 10 23 cm -2 ) of the gas column density from 
molecular absorption measurements in the radio (Wiklind & Combes 1995; Gerin et al. 1997; 
Carilli, Rupen & Yanny 1993). We considered only the two brightest images in each system, 
and for Q 2237+0305, where several extinction measurements are available, we used the value 
from Falco et al. (1999) to be consistent with other lenses. The values for Q 2237+0305 are 
consistent with those measured by Agol, Jones & Blaes (2000). For estimating the average 
dust-to-gas ratio we must exclude HE 1104—1805 (see §3.2). If we fit the four remaining 
lenses for the mean gas-to-dust ratio including the uncertainties in both the absorption and 
the extinction, as well as allowing for additional intrinsic scatter in the ratio beyond the 
measurement errors, we find a mean ratio of (1.4 ± 0.5) x 10 -22 mag cm 2 atoms -1 and 
estimate that there is an intrinsic scatter in the ratio of ~ 40% beyond our measurement 
errors. If we allow for no intrinsic scatter we find a ratio of (1.5 ± 0.3) x 10 -22 mag cm 2 
atoms -1 . These estimates are consistent with the average Galactic value of 1.7 x 10 -22 
mag cm 2 atoms -1 but show a modestly larger intrinsic scatter than the locally observed 30% 
(Bohlin, Savage & Drake 1978). The estimates of the intrinsic scatter is consistent with local 
estimates given our small sample size. Figure 4 compares the results for the lens galaxies to 
the local stellar sample from Bohlin, Savage & Drake (1978). 

The lens galaxies of our small sample consist of three spiral and one elliptical galaxies, 
although the quasar in Q 2237+0305 is lensed by the central bulge of the spiral galaxy 
(Table 2). Considering the environmental and redshift differences, it is interesting that we 
obtained dust-to-gas ratios consistent with those of our Galaxy. It suggests that the dust- 
to-gas ratio of the ISM is similar in all normal galaxies. It may not, however, be universal, 
since Maiolino et al. (2001) found that the circumnuclear regions of local AGN may have a 
very low dust-to-gas ratio (a factor of ~ 3-100 lower than Galactic), possibly because of grain 
formation and destruction processes peculiar to the local environment created by the AGN 
(Maiolino et al. 2001; Maiolino, Marconi & Oliva 2001). Recently, models of evolution of 
dust-to-gas ratios are discussed in several papers (e.g., Dwek 1998; Edmunds 2001; Inoue 
2003). However, we are unable to compare the details of the models with our small sample 
of data. 

It is relatively straight forward to improve on these results by systematically obtain- 
ing resolved X-ray spectra of lenses showing significant differential extinctions. Like most 
X-ray spectral studies, the required integration times are fairly long (~ 50 ksec). The dif- 
ferential absorption measurement requires resolved spectra from individual lensed images. 
This means that XMM-Newton cannot be used despite its higher sensitivity - the high 
spatial resolution of Chandra is required. For systems with modest differential extinctions 
(AE(B — V) < 0.1 mag) there are probably significant systematic uncertainties created by 
chromatic microlensing where color differences created by microlensing are misinterpreted 
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as extinction (e.g., Wucknitz et al. 2003). The X-ray emitting region is probably not as 
spatially extended as the optical accretion disk, so that the entire X-ray continuum is mag- 
nified by roughly a constant at a given epoch, unless the observations occur during a high 
micro lensing magnification event. Thus, the column density estimates obtained from the 
absorption in the soft X-ray spectra should be insensitive to microlensing effects. Focusing 
on lenses with significant differential extinction will also minimize the possibility that sig- 
nificant differences in the properties of the ISM between lensed images will skew the results 
either because the extinction laws are very different (e.g., Munoz et al. 2004; McGough et al. 
2005) or because the dust-to-gas ratios are very different. With a large enough sample, it 
would be possible to study the evolution of the ISM with redshift. 
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